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ABSTRACT: Colliding beams of 70% polarized protons at up to ^=500 GeV, with high 
luminosity, L=2xl0^^ cm~'^sec~^, will represent a new and unique laboratory for studying the 
proton. RHIC-Spin will be the first polarized-proton collider and will be capable of copious 
production of jets, directly produced photons, and W and Z bosons. Features will include 
direct and precise measurements of the polarization of the gluons and of u, d, u, and d quarks 
in a polarized proton. Parity violation searches for physics beyond the standard model will be 
competitive with unpolarized searches at the Fermilab Tevatron. Transverse spin will explore 
transversity for the first time, as well as quark-gluon correlations in the proton. Spin dependence 
of the total cross section and in the Coulomb nuclear interference region will be measured at 
collider energies for the first time. These qualitatively new measurements can be expected to 
deepen our understanding of the structure of matter and of the strong interaction. 
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1 INTRODUCTION 



Spin is a powerful and elegant tool in physics. One of the most exciting aspects of 
physics is a search for the unexpected, the nonintuitive, in nature. Intrinsic spin 
itself violates our intuition, in that an elementary particle such as an electron can 
both be pointlike and have a perpetual angular momentum. We find at this time 
an apparent violation of our intuition in the proton. We understand the proton as 
being composed of quarks, gluons, and antiquarks, and we expect the proton spin 
to be carried dominantly by its three valence quarks. Instead, through the 1980s 
and 1990s, deep inelastic scattering (DIS) experiments of polarized electrons and 
muons from polarized nucleons have shown that on average only about 1/4 to 
1/3 of the proton spin is carried by the quarks and antiquarks in the proton [|| . 
Therefore, the spin of the proton appears to be mainly carried by the gluons and 
orbital angular momentum! This surprising and counterintuitive result indicates 
that the proton, and particularly its spin structure, is much more interesting than 
we had thought. 

Spin can be used as an elegant tool to search for the unexpected. If an exper- 
iment is found to depend on the spin direction, it can violate a deep expectation 
that physics should be symmetric with respect to that axis. An example is mirror 
symmetry, that physics should not depend on left- or right-handedness. The vi- 
olation of parity by the weak interaction was the surprise that led to the present 
electroweak model with the purely left-handed charged weak vector bosons W^. 
At the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Labora- 
tory, the and W~ will be produced by colliding beams of protons spinning 
alternately left- and right-handed. The expected maximum violation of parity 
will allow unique and precise measurements of the spin direction of the quarks 
and antiquarks in the proton that form the W bosons, identified by quark fla- 
vor, u, u, d, and d. A dependence on handedness in the production of jets at 
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RHIC beyond the contribution from W and Z would directly signal new physics, 
possibly coming from quark substructure at a scale above the weak scale. 

Physics is also a search for unexpected order in nature. Large spin effects 
necessarily imply coherence and order. If the gluons in a proton are found to 
be dominantly spinning in the same direction, as discussed widely in the context 
of the smallness of the quark spin contribution (reviewed in Q), there would 
need to be a simple underlying physical mechanism that creates this order. At 
RHIC, dedicated experiments will measure the direction of the gluon spin in the 
proton for the first time — an exciting prospect, since there are hints that the 
gluon polarization may be substantial. 

The RHIC at Brookhaven has begun a program of colliding beams of gold ions 
at 100 GeV per nucleon in the spring of 2000. The following year, the first physics 
run colliding beams of polarized protons is expected. RHIC-Spin will be the first 
polarized proton-proton collider. It will reach an energy and luminosity at which 
the collisions can clearly be interpreted as collisions of polarized quarks and 
gluons, and it will be capable of copious production of jets and directly produced 
photons, as well as W and Z bosons. Quantum chromodynamics (QCD) makes 
definite predictions for the hard spin interactions of quarks and gluons, which 
implies that RHIC will enable us to test a sector of QCD that so far has been 
little explored. The polarized quark and gluon probes at RHIC complement the 
beautiful work done using polarized lepton probes to study proton spin structure. 
These strong interaction probes will be sensitive to the gluon polarization in jet 
and direct photon production and will allow quark spin-flavor separation in 
production. RHIC-Spin will also represent the highest energy for proton-proton 
collisions at accelerators, and unpolarized production will be used to precisely 
measure the flavor asymmetry of the antiquark sea. 

At the Polarized Collider Workshop at Penn State University in 1990 j^, the 
exploration of the spin of the proton was a major focus for the physics of polarized 
proton collisions at RHIC. The RHIC Spin Collaboration was formed the follow- 
ing year, consisting of experimenters, theorists, and accelerator physicists |^]. 
Since 1993, the two large heavy ion detectors at RHIC, Star and Phenix, have 
considered spin as a major program and include additional apparatus specifically 
for spin physics. In addition, the pp2pp experiment at RHIC, studying small- 
angle elastic scattering, will also feature spin. The present article presents the 
anticipated physics of the RHIC spin program as developed by the RHIC Spin 
Collaboration and by the Star, Phenix, and pp2pp Collaborations. 

The RHIC spin accelerator complex is illustrated in Figure 1. An intense 
polarized H~ source feeds a chain of accelerators. Individual bunches of 2 x 10^^ 
protons with 70% polarization are transferred from the Alternating Gradient 
Synchrotron (AGS) to the RHIC rings at 22 GeV. This is repeated 120 times for 
each ring at RHIC. The polarized protons are then accelerated to up to 250 GeV 
in each ring for collisions at each of 6 intersection regions. With a f3* = 1-meter 
focus at Star and Phenix, luminosity will be £ = 2 x lO'^^ cm~^s~^, for the 
highest RHIC energy of -^s = 500 GeV. Experimental sensitivities given in this 
article are based on 800 pb^^ for ^ = 500 GeV and 320 pb"^ for ^ = 200 GeV. 
This corresponds to runs of 4 x 10^ s at full luminosity, about four months of 
running with 40% efficiency, at each energy. We expect the data to be collected 
over three to four years, since RHIC is shared between heavy-ion and polarized- 
proton collisions. The expected sensitivities will be excellent due to the high 
luminosity for proton-proton collisions. For comparison, we note that the pp 
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Tevatron at Fermilab has run for a total of ~ 130 pb^i as of 1999. 

It is difficult to maintain the proton polarization through acceleration because 
of its large anomalous magnetic moment: the proton spin readily responds to 
focusing and error magnetic fields in the rings, and spin resonances are encoun- 
tered frequently, for example at every 500 MeV of acceleration in the AGS. The 
methods that are used to avoid depolarization in acceleration are very elegant, 
and the acceleration of polarized protons to 250 GeV will be breaking new ground 
in accelerator physics. The key device is a string of dipole magnets that rotate 
the proton spin 180° around a selected axis in the horizontal plane each time 
the beam passes H]. Each two passes in effect cancel the cumulative tilt of the 
spin resulting from horizontal magnetic fields, thus eliminating the major spin 
resonances at RHIC. There will be four "Siberian Snakes" at RHIC, two in each 
ring. The name refers to the home institution of the inventors (Novosibirsk) and 
to the motion of the beam passing through. In this article, we do not discuss 
the accelerator physics work leading to the RHIC spin plan but, as for any 
spin experiment, past or future, there is a very tight, necessary, and refreshing 
coupling between the polarization technology and the physics. 

For two Siberian Snakes in each ring, the stable spin direction in RHIC will 
be vertical. Therefore, transverse spin physics will be available to all the exper- 
iments. For Star and Phenix, special strings of dipole magnets will be used 
to rotate the spin to longitudinal at their intersection regions. Longitudinal spin 
is necessary to study gluon polarization and parity- violating physics. A recent 
plan is to initially use one Siberian Snake in each ring, which allows the con- 
struction and installation of the Snakes and Rotators to be staged. With a single 
Snake in a ring, the stable spin direction is in the horizontal plane. If the beam 
is inserted into RHIC, and the Snake is then turned on adiabatically, the spin 
will follow from vertical to horizontal. At energies roughly 2 GeV apart, it will 
be possible to have longitudinal polarization at all six intersection regions, up to 
a beam energy of 100 GeV. One Snake is already installed in RHIC at this time, 
and a second Snake will be completed in summer 2000. Therefore, the RHIC-Spin 
program will be ready for its commissioning in summer 2000 and ready for the 
first spin physics run with longitudinal polarization at y/s = 200 GeV in 2001. 
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Polarized Proton Collisions at BNL 




Figure 1: Schematic layout of the RHIC accelerator complex. Only relevant devices for 
polarized pp collisions are shown. 
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2 PREREQUISITES FOR SPIN PHYSICS AT RHIC 

2.1 Theoretical Concepts and Tools 

2.1.1 Partons in High-Energy Scattering: Factorization 

Polarized pp collisions at RHIC will take place at center-of-mass energies of 
^/s =200-500 GeV. Except for polarization, we have a typical collider physics 
situation, similar to that at CERN's SppS or the Tevatron at Fermilab. One 
therefore expects that parton model concepts, augmented by the predictive power 
of perturbative QCD, will play a crucial role in describing much of the interesting 
spin physics to be studied at RHIC, if the reaction under consideration involves 
a hard probe, for instance a photon produced at transverse momentum (pt) of a 
few GeV or more. 

The QCD-improved parton model has been successfully applied to many high- 
energy processes involving hadrons in the initial or final state. In this framework, 
a cross section is written in a factorized form as a convolution of appropriate 
parton densities and/or fragmentation functions with a partonic subprocess cross 
section. The predictive power of perturbative QCD follows from the universality 
of the distribution functions: Once extracted from the data in one process, they 
can be used to make definite predictions for any other. As an example, let us 
consider the production of a pion with large pT in a collision of unpolarized 
protons, that is, pp irX. The process is depicted in Figure |2|. In the parton 




Figure 2: Production of a large-pr pion in a hard pp collision. 

model framework, in the context of QCD perturbation theory, one writes the 
cross section convolution, 



Jdxidx2dzff{xi,i/)f^{x2,n'^ 



/l,/2,/ 

X {xiPi,X2P2,P7t/z,h)DJ{z,^ ), (1) 

where pi and p2 are the incident proton momenta. Here, V stands for any ap- 
propriate set of the kinematic variables of the reaction. Furthermore, ff (x, p,'^) is 
introduced as the probability density for finding a parton of type fi in the proton, 
which has taken fraction x of the proton's momentum. Likewise, -CJ (-2, /^^) is the 
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probability density for finding a pion with momentum fraction z in the parton 
/. The (T-''i'^2^/^ gj-Q ^j^g underlying hard-scattering cross sections for initial 
partons /i and /2 producing a final-state parton / plus unobserved X'. 

The functions and DJ introduced in Equation |^ express intrinsic proper- 
ties of the proton and of the hadronization mechanism, respectively. Therefore, 
they are sensitive to non-perturbative physics and cannot be calculated from first 
principles in QCD at present. In contrast to this, for a sufficiently hard process, 
it will make sense to calculate the subprocess cross sections af^S'^^f^ as per- 
turbation series in the strong coupling Ug- The separation of short-distance and 
long-distance phenomena as embodied in Equation || necessarily implies the in- 
troduction of an unphysical mass scale ^, the factorization scale. The presence 
of jjL arises in practice when computing higher-order corrections to the af^S'^^f^ . 
Here, one encounters singularities resulting from configurations in which one of 
the incoming (massless) partons collinearly emits another parton. In the same 
way, such "coUinear" singularities (or "mass" singularities) occur in the final 
state from collinear processes involving parton /. Regularization of the mass 
singularities always introduces an extra mass scale M to the problem; the cross 
section depends on it through powers of "large" logarithms of the type ln(pr/M). 
The collinear-singular logarithms are separated off at the factorization scale /i, 
to be of the order of the hard scale pT characterizing the hard interaction, and 
are absorbed ( "factorized" ) into the "bare" parton densities (or fragmentation 
functions). This procedure is of use only if it is universal in the sense that the 
mass singularities absorbed into the parton densities are the same for all pro- 
cesses involving a given initial parton. Proof of this property is the subject of 
factorization theorems ^ and is necessary for the parton model to be valid in 
the presence of QCD interactions. 

In summary, the QCD-improved parton-model picture as used for Equation 1 
consists of perturbatively calculable partonic hard-scattering cross sections and of 
scale-dependent parton densities and fragmentation functions that are universal 
in the sense that once they are measured in one process, they can be used to 
make predictions for any other hard process. It is important to point out that the 
parton densities and fragmentation functions are never entirely nonperturbative: 
Their dependence on the factorization scale is calculable perturbatively, once 
the densities are known at some initial scale /io- This has to be so, since the 
/i-dependence of the (jf^f'^^f-^ is calculable and the prediction of a physical 
quantity, such as the hadronic cross section a^P^'"^ ^ has to be independent of 
IJL to the order of perturbation theory considered. The tool to calculate the 
dependence of the and D'j on the "resolution scale" is the set of evolution 
equations [10|. 



2.1.2 Spin-Dependent Parton Densities and Cross Sections 

So far we have disregarded the spin information contained in parton distributions 
and fragmentation functions. If a hard-scattering process with incoming protons 
having definite spin orientation is studied, as at RHIC, one expects it to give 
information on the spin distributions of quarks and gluons in a polarized proton. 
The possible parton distribution functions are summarized in Table |^. A 
similar table could be presented for polarized fragmentation functions [^]: The 
observation of the polarization of a final-state hadron should give information on 
the polarization of the parton fragmenting into that hadron. 
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Table 1: Compilation of quark and gluon parton densities including spin depen- 
dence. The ubiquitous argument {x,fi'^) of the densities has been suppressed. 
For brevity, a column for antiquarks (q) was omitted, which would have an iden- 
tical structure to that of the quark column. Labels -|-, — denote helicities, and 
t, i transverse polarizations. Superscripts refer to partons and subscripts to the 
parent hadron. 



Polarization 


Quarks 


Gluons 


unpolarized 


q = ql + q~=qj+q^ 


g = 9l+9+ 


long, polarized 




^9 = 9+- 9+ 


transversity 


oq = q\ - q\ 





Within roughly the past decade, beautiful data have become available that 
are sensitive to the "longitudinally" polarized ( "helicity- weighted" ) parton den- 
sities of the nucleon. The tool to obtain such information has been deep-inelastic 
scattering (DIS) of longitudinally polarized leptons and nucleons. The spin asym- 
metry measured in such reactions gives information on the probability of finding 
a certain parton type {f = u,u,d,d, . . . , g) with positive helicity in a nucleon of 
positive helicity, minus the probability for finding it with negative helicity (see Ta- 
ble H). These densities are denoted as A/(x,/x^). The Appendix provides a brief 
discussion of the implications of present polarized DIS data on our knowledge 
about the A/. Within a parton-model concept, the integrals of the A/(x,/x^) 
over all momentum Bjorken-x ("first moments"), multiplied by the spin of the 
parton /, will by definition give the amount of the proton's spin carried by species 
/, appearing in the proton-spin sum rule: 



dx 



1 



Y,{^q + ^q){x.^^^) + ^9{x,fi^ 



where L is the orbital angular momentum of quarks and gluons in the proton |13|. 

The longitudinally polarized parton distributions in Table |l] can be separated 
from the unpolarized ones if suitable differences of cross sections for various lon- 
gitudinal spin settings of the initial hadrons are taken [14|: 

tX r pp—*ttX J pp^ttX j pp—*nX 



dAaPP- 



tX 



dV 



da 



da 



dV 



dV 



dV 



+ 



da- 



dV 



^ dxidx2dz Af^{xi,fi'^) Af^{x2,n'^) 



dV 



-{xi,pi,X2,P2,Ptt/z,ij) D]{z,p ) , 



where 

dV 



da_ 



;,/l/2-/X' ^^hf2^fX' 



+- 



da{ 



dV 



dV 



dV 



■ + ■ 



dai 



dV 



(3) 



• (4) 
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Here and in Equation ^ subscripts denote the helicities of the incoming parti- 
cles, i.e. of the protons in Equation ^ and of partons /i, /2 in Equation HI Thus, 
the "longitudinally polarized" cross section dAa^P^'"-^ /dP depends onlyQ on the 
parton densities for longitudinal polarization and on the (calculable) "longitu- 
dinally polarized" subprocess cross sections dAa-^^-^^^/^ /dV. A measurement 
of dAaPP^^'^/dV therefore gives access to the A/. Adding, on the other hand, 
all terms in the first line of Equation ^, one simply returns to the unpolarized 
cross section in Equation ^, with its unpolarized densities / and the unpolarized 
subprocess cross sections d(jf^f'^~*f^ /dV, corresponding also to taking the sum 
of the terms in Equation |^. 

Notice that we have taken both initial protons to be polarized in Equation 
If only one is polarized, we can still define a singly polarized cross section by 
^^PP^TT i^j) _ ^^pp^TT jfjij)^ where the subscript refers to the polarized proton's 
helicity. However, this combination can be nonzero only if parity is violated in 
the hard process [0]. If so, the single-spin cross section will depend on products 
of parton densities A/i and /2, representing the polarized and the unpolarized 
proton, respectively. 

With two transversely polarized beams, one will take the first line of Equation^ 
for transverse polarizations rather than helicities. The result will be a polarized 
cross section depending on transversely polarized subprocess cross sections and, 
for each proton, on the differences of distributions of quarks (or antiquarks) 
with transverse spin aligned and anti-aligned with the transverse proton spin. 
The latter quantities are the "transversity" distributions 11, 16, 17] and are 
denoted 6f{x, fi^) (see Table |^).[| Note that in the case of transverse polarization a 
cos(2(/)) dependence of the cross section on the azimuthal angle (j) of the observed 
final-state particle arises |15, 11, |l^, 17], since an extra axis is defined by the 
transverse spin. We also mention that transverse single-spin cross sections, such 
as da^^^ /dV - daf^"" /dV, are allowed to be nonzero in QCD but vanish in 
the simple parton- model picture presented so far ]|l^, |l9| (see Section ^ ). 

Extension to polarization in the final state is also possible. If the observed 
particle in Equation ^ were, say, a A-hyperon instead of the (spinless) pion, one 
could consider the first line of Equation ^ for the helicities of one of the incoming 
protons (the other proton is assumed to be unpolarized, for simplicity) and of the 



A. In this way one obtains a "helicity transfer" cross section |14] that depends on 
the distribution of parton /2 for the unpolarized proton, on A/i for the polarized 
proton, on polarized fragmentation functions ADj (defined in analogy with A/), 
and on helicity-transfer subprocess cross sections. 

For spin experiments, the most important quantity in practice is not the po- 
larized cross section itself, but the spin asymmetry, which is given by the ratio of 
the polarized over the unpolarized cross section. For our example above, it reads 

dAaPP^-^/dV 
^LL daPP^^^/dV ■ 

For the asymmetry, one often uses subscripts to denote the type of polarization 
(L=longitudinal, T=transverse) of the initial particles. As follows from Equa- 

^In addition, there is dependence on the pion fragmentation functions DJ . 
^One frequently also finds the notation ATf{x,n^) or h({x,^^) in the Uterature. 
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tion 1^, the resulting spin asymmetry will possess the generic structure 



A 



LL 







xDj 




X 







where d' 



hf2-*fx' 

LL 



(6) 



d^a^'^f'^^f'^' /da^^f^^f-^' is the spin asymmetry for the sub- 



process /1/2 fx' , often also referred to as the analyzing power of the reaction 
considered. The lowest-order analyzing powers for many reactions interesting at 
RHIC are depicted in Figure |3|. 
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Figure 3: Lowest-order analyzing powers for various reactions relevant for RHIC, as 
functions of the partonic center-of-mass system (cms) scattering angle Left: 
longitudinal polarization, right: transverse polarization (a factor cos(2(/>) has been taken 
out, where is the azimuthal angle of one produced particle). 



2.2 Detection 

2.2.1 Asymmetries and Errors 

Asymmetries in a collider experiment can be defined (and measured!) for a single 
polarized beam or for both beams polarized, with longitudinally polarized beams, 
transversely polarized beams, or with a combination of these. Additionally, one 
can study a combination of beam spin state and final-state angular dependence. 
For longitudinal polarization for both beams, the asymmetry All is defined as 

^ ^ + (7—) - (o-+- + <y-+) 

(cr++ + cr__) + (o-+_ + cr_+) 



(7) 
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Here, (T_| represents a cross section for producing a specified final state witli the 

initial proton helicities (+) and (— ). However, the proton beams are not in pure 
helicity states. We expect that the beams will be about 70% polarized, meaning 
that 

where refers to the number of protons in the beam with (+) helicity. There- 
fore, collisions with two bunches of protons, with for example +0.7 polarization 
for one bunch and —0.7 polarization for the other bunch, will include collisions of 

all four helicity combinations, (++), (H — ), ( — h), and ( ). The experimental 

asymmetry is defined as follows: 

A,, ' ^, (iv;+ + ivi-)-(iv;- + ivi+) .g. 
P1P2 (iv;+ + ivi_) + (iv;_ + ivi+)' ^ ^ 

where N'_^_ represents the observed number of events when the beams were po- 
larized (+) for beam 1 and (— ) for beam 2, and normalized by the luminosity 
for the crossing. Here, it is necessary to know only the relative luminosity for 

the (++) and ( ) collisions versus the (H — ) and ( — h) collisions. The beam 

polarizations are Pi and i-*2- Algebra can confirm that Equation ^ is equivalent 
to Equation |^. 

Similarly, we can define the parity-violating asymmetry for one beam polarized 
longitudinally, 

. 1 N'-N'_ 
Al = — , Al = X -. (10) 

The parity- violating asymmetry was defined in 1958 to be positive for left-handed 
production Observed parity- violating asymmetries are therefore typically 

positive, due to the left-handed weak interaction. 

For transverse spin, one- and two-spin asymmetries are defined in analogy with 
the longitudinal asymmetries above, referred to as and Att- In this case, the 
directions (+) and (— ) are transverse spin directions of the beam protons, not 
the helicities. The transverse-spin asymmetries depend on the production angle, 
9, and on the azimuthal angle of the scattering, cj), as well as other variables. The 
azimuthal dependence for scattering two spin-1/2 particles is 

Att oc cos(2</)) and An oc cos((/>). (11) 

(^=0 is defined for scattering in the plane perpendicular to the polarization di- 
rection. Typically the beam is polarized vertically, with (+) polarization up, and 
positive A]\j implies more scattering to the left than to the right of the beam 
direction. The notation Ann is also used for a transverse two-spin asymmetry, 
where refers to beam polarization normal to the scattering plane. A subscript 
S traditionally designates beam polarization in the transverse direction in the 
scattering plane. 

Prom Equation ^ or Equation |l^ we need to know the beam polarization(s), 
count the number of signal events for each combination of beam spin directions, 
and monitor the relative luminosity for the crossings with these combinations of 
beam spin directions. The statistical error of the measurement is 
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Here N is the total number of events observed, and it is assumed that the statis- 
tical errors on the relative luminosities and on the beam polarization are small. 
For the single spin asymmetry, 

For small to moderate asymmetries, 

AAll = ±l/iPiP2) X ^ and AAl = ±1/P x (14) 

V jV y Jy 

Since we expect P = Pi = P2 = 0.7, 10^ events would give an error of AAll = 
±0.02 for the double spin asymmetry, or AAl = ±0.014 for the parity- violating 
asymmetry. 

In principle, asymmetry measurements are very straightforward. As long as 
the detector acceptance remains stable with time between reversals of the beam 
spin states, the measurement will be stable and the errors will be largely sta- 
tistical. However, when reversals of the beam polarization are spread apart in 
time, and/or the beam conditions for opposite spin states differ, acceptance can 
change and false asymmetries develop. At RHIC the bunches, 120 in each ring, 
are prepared independently at the source, so that the bunches can alternate po- 
larization sign, 106 ns apart, as shown in Figure Note that one ring with 
alternate bunches and the other ring with alternating pairs of bunches create 
the four spin combinations, (±±), (H — ), ( — h), and ( ). Therefore, the con- 
cern of time-dependent acceptance and beam location variations for opposite sign 
beams should be negligible at RHIC, and asymmetry measurement errors should 
be mainly statistical, even for small asymmetries. 




Figure 4: Bunch filling pattern with respect to the spin states of polarized protons. 

What systematic errors do we expect at RHIC? There are two classes of sys- 
tematic errors: false asymmetries and scale errors. If the relative luminosities 
for the bunch spin combinations are incorrectly measured through, for exam- 
ple, a saturation effect in the luminosity monitor, which couples to variations in 
beam intensity for the bunch spin combinations, the numerator of Equation |^ or 
Equation ^ will be nonzero from the incorrect normalization, creating a false 
asymmetry. If the beam polarization is incorrect, no false asymmetry is created, 
but the scale of the resulting asymmetry is changed. 

Each experiment will measure the relative luminosities for each crossing. The 
luminosity monitors must be independent of beam polarization, and statistical 
errors on the relative luminosity measurements need to be very small to match 
the statistical sensitivity available for high-statistics measurements, such as jet 
production. 
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Relative luminosity needs to be known to the 10~^ level for some asymmetry 
measurements. This job appears daunting, but the time dependence of the accep- 
tance (efficiency is included with acceptance in this discussion) of the luminosity 
monitor needs to be stable only over roughly one turn of RHIC, or 13 /.is. 

2.2.2 Polar imetry 

Polarization is measured by using a scattering process with known analyzing 
power. Knowledge of the analyzing power for different processes can come from 
theoretical calculations, for example for QED processes, and from experimental 
measurements using a beam or target with known polarization. Polarimetry 
at RHIC will be based on elastic proton-proton and elastic proton-carbon 
scattering in the Coulomb nuclear interference (CNI) region. The analyzing power 
there is largely calculable; it is expected to be small but significant. It can be 
determined to excellent precision using a polarized proton target in RHIC, and 
the rates for CNI scattering are very high. 

Sensitivity to the proton spin is from scattering the Coulomb field of an unpo- 
larized particle (proton or carbon) from the magnetic moment of the polarized 
proton. This method uses the dominance of the interference of the one-photon ex- 
change helicity-flip electromagnetic amplitude, proportional to the proton anoma- 
lous magnetic moment, with the non-fiip strong hadronic amplitude, which is 



determined by the pp or pC total cross section fitot [Hi H^i ^3 HI, However, 
there can also be a hadronic spin-flip term, which is not presently calculable. 
(This possibility is discussed in more detail in Section 7.) Therefore, significant 
sensitivity to the proton spin is predicted over the entire RHIC energy range from 
the electromagnetic term, but the absolute sensitivity is limited to ±15% [p5|] . 
For this reason, a polarized hydrogen gas jet target will be installed at RHIC. 
The polarization of the jet target can be measured to ±3% so that the analyzing 
power in the CNI region can be measured precisely, and this analyzing power will 
then be used to determine the beam polarization at RHIC precisely. 

Existing polarized hydrogen gas jet targets are thin, so that the determination 
of the beam polarization using the jet target will take hours. For this reason, 
RHIC will also use carbon ribbon targets and use proton-carbon CNI scattering 
to monitor the beam polarization frequently. 

Absolute beam polarization is expected to be known to ±5%. The systematic 
scale uncertainty of the asymmetry measurements will be of the order of ±5% 



for single spin measurements such as Equation 10 and ±10% for two spin mea- 
surements such as Equation |9|. By scale uncertainty we mean that in forming 
the ratio of the error in the asymmetry measurement, A All in Equation 0, 
to the measurement itself. An in Equation ^, the polarization normalization 
divides out. Therefore, the polarization uncertainty applies to the scale of the 
measurement and not to the statistical significance of the measurement. 

2.2.3 RHIC Detectors 

This article emphasizes the physics that will be probed at RHIC-Spin. There 
are six collision points at RHIC, as shown in Figure |l|, and two are used for the 
two large detectors, Phenix [|^] and Star |28|. These detectors are quite com- 



plementary: Star emphasizes large coverage with tracking, and the strengths of 
Phenix are in fine-grained calorimetry for photons and electrons and in "forward" 
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muon detectors. Sensitivities for the spin measurements at RHIC are based on 
these detectors. Although one could discuss the sensitivity for a 47r-acceptance 
fine-grained detector, such a detector does not exist. And we note that, for 
example, the Phenix electromagnetic calorimeter (EMCal) has 100 times finer 
granularity than previous collider detectors. The pp2pp ||2^ and Brahms 
detectors share one collision point, and the Phobos detector is located at 
another crossing. Star and Phenix will measure gluon and quark polarizations 
with hard scattering. The pp2pp experiment will measure spin dependence in 
small-angle elastic scattering; Brahms and Phobos will measure transverse spin 
asymmetries. 




Figure 5: The Phenix detector system. 

The Phenix detector, shown in Figure ||, has two central arms at 90° to the 
beams with fine-grained EMCal towers, Ar] x A(/> = 0.01 x 0.01. The minimum 
opening angle for vr*^ 77 corresponds to one tower for a 30-GeV vr*^. This is 
important to separate directly produced photons, a probe of gluon polarization, 
from background from decay. Resolution is excellent, with AE/E = ±3% 
at 10 GeV. The two central arms each cover 90° in azimuth, left and right. 
Pseudorapidity acceptance is |?7| < 0.35. The vertex detector, central tracker, 
ring-imaging Cherenkov detector (RICH), and time expansion chamber (TEC) 
are also shown. The central magnetic field, provided by two Helmholtz coils, is 0.8 
Tesla meters, integrated radially. The tracking pT resolution is Apx/px = ±2.5% 
at 10 GeV/c for the east arm, which includes the TEC, and ±5% in the west 
arm without a TEC. Triggering in the central arms, for selection of high-pj^ direct 
photons, electrons, and charged pions, will be based on overlapping tower clusters 
in the EMCal, combined with RICH information. Studies indicate a sufficiently 
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clean and efficient electron trigger to allow pT > ^ GeV/c or so. Such a low- 
momentum electron trigger is attractive to obtain charm quark events. 



STAR Detector 

Silicon Vertex 




Figure 6: Cut view of Star detector system. 

The Phenix muon arms surround the beams, covering A(j) = 2tt and 1.2 < 
\r]\ < 2.4. The arms include a muon identifier (MuID) with five iron-detector 
layers, as well as three tracking stations. The muon arm magnets produce a 
radial field, ranging from 0.2 to 0.75 Tesla meters, integrated along the beam 
direction. Longitudinal momentum resolution is about ±2% at 10 GeV/c. A 
4-GeV/c muon penetrates to the fifth MuID layer. 

Phenix will emphasize muon measurements for W — > fiv, Drell-Yan lepton 
pairs, J/V') and heavy fiavors. Central arms will measure 7, jet fragmentation to 
TT*^'^, and W eu, as well as heavy fiavors (single lepton, and e with p.), with 
small acceptance and high granularity. 

The Star detector is shown schematically in Figure ^. A barrel time projection 
chamber (TPC) covers \ri\ < 1.0 and Acp = 2ti. This is surrounded by an EMCal 
with towers Ar/ x A.(f) = 0.05 x 0.05. The EMCal has a shower maximum detec- 
tor at a depth of five radiation lengths with projective readout wire chambers 
reading longitudinally and azimuthally, each with 1-cm spacing. Studies show an 
effective separation of single photons from merged photons from 7r° decay out to 
Pt = 25 GeV/c. Energy resolution is excellent, with IS.E/E = it 5% at 10 GeV. 
Additional barrel detection includes a silicon drift vertex tracker around the colli- 
sion point and an array of trigger counters outside the TPC. The central solenoid 
field is 1.0 Tesla meters, integrated radially. The Star px resolution is ±3% at 
Pt = 10 GeV/c. Star is also building one endcap calorimeter to cover 1 < 77 < 2 
for photons and electrons and to expand the jet cone coverage. 
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Triggering at Star will be based on the trigger counters and EMCal, which are 
fast detectors. A major issue to resolve is the long memory of the TPC, which will 
include on the order of 800 out-of-time tracks from the 40-/Lts drift time, at full 
luminosity with a 10-MHz collision rate. Star studies have shown that the drift 
of the out-of-time tracks cause them to point significantly away from the collision 
point. This drift will be used to remove the unwanted tracks, and this must be 
done before writing to tape. Studies have also shown good jet reconstruction after 
the tracks are removed. Jets will be reconstructed at Star with a combination 
of EMCal and tracking, with no hadronic calorimetry. Simulations show a full 
width at half maximum of 30% for the Apjet/pjet distribution, limited by the 
hadronization dynamics of final-state partons. A cone size of i? = ^ Arj"^ + A(j)'^ = 
0.7 was used. Star will measure jets, 'j+jet, and W — > ev, with wide acceptance 
and reconstruction of the parton kinematics. 

The PP2pp experiment is designed to study small-angle proton-proton elastic 
scattering, from —t = 0.0005 to —t = 1.5 (GcV/c)^. Silicon-strip detectors will 
be placed in Roman pots at two locations along each beam to measure scattering 
to very small angles. The experiment will also use a polarized hydrogen jet target 
with silicon recoil detectors to cover lower center-of-mass energy and will measure 
the absolute polarization of the RHIC beams. 

The Brahms detector has two movable spectrometers (2.3° < < 30° and 
30° < < 95°) with superb particle identification. The spectrometer covers up 
to 30 GeV/c with Ap/p = ±0.1% and it will provide unique measurements of 
single transverse-spin asymmetries in the forward, thus high-xj?, region. 

Phobos is a table-top-sized detector that uses silicon-strip detectors to cover 
a large solid angle. Two spectrometers comprise 15 planes of silicon pad detectors 
each, with 7 planes in a 2-Tesla magnetic field. Its wide geometrical acceptance 
and momentum resolution is suitable for pair or multiparticle final states in spin 
physics, such as p° — tt+tt". 
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3 MEASURING Ag AT RHIC 

Measurement of the gluon polarization in a polarized proton is a major emphasis 
and strength of RHIC-Spin. By virtue of the spin sum rule (|^), a large Ag is 
an exciting possible implication |^ of the measured |j^ smallness of the quark 
and antiquark contribution to the proton spin. A large gluon polarization would 
imply unexpected dynamics in the proton's spin structure. Because of this special 
importance of Ag, and since it is left virtually unconstrained by the inclusive-DIS 
experiments performed so far (see Appendix), several experiments focus on its 
measurement. A fixed-target DIS experiment, Hermes, measures the process 



e{^)p — > h^h~X |32|, where h = 7r,K, which is in principle sensitive to the gluon 
polarization. However, the transverse momenta are low, making interpretation 
in a hard-scattering formalism difficult. The DIS experiment Compass (see e.g. 
|l33|| ) will measure the same reaction at higher energies, as well as heavy-fiavor 
production, to access gluon polarization. Scaling violations and the reaction 
e{'y)p — > jet(s)X will constrain at HERA, if the proton ring is polarized [^]. 
At RHIC, the gluon polarization will be measured directly, precisely, and over a 
large range of gluon momentum fraction, with large momentum transfer ensuring 
the applicability of perturbative QCD to describe the scattering, and with several 
independent processes. The RHIC probes, shown in Figure ^, are as follows: 

• High-pT ( "prompt" ) photon production pp 

• Jet production, pp jet(s)X 

• Heavy-flavor production, pp ccX, bbX 



(a) (b) (c) 




Figure 7: Selected lowest-order Feynman diagrams for elementary processes 
with gluons in the initial state in pp collisions: (a) quark-gluon Compton process 
for prompt-photon production, (6) gluon-gluon and gluon-quark scattering for jet 
production, and (c) gluon-gluon fusion producing a heavy quark pair. 



3.1 Prompt- Photon Production 



Prompt-photon production, pp,pp,pN jX [35], has been the classical tool for 
determining the unpolarized gluon density at intermediate and large x. At leading 
order, a photon in the final state is produced in the reactions qg — > (Figure |^a) 
and qq ^g. Proton-proton, as opposed to proton-antiproton, scattering favors 
the quark-gluon Compton process, since the proton's antiquark densities are much 
smaller than the quark ones. The analyzing power for direct photon production is 
large (Figure S). Photons produced in this way through partonic hard scattering 
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show a distinct signal at colliders, that of an isolated single photon without jet 
debris nearby. The production of photons with polarized beams at RHIC is 
therefore a very promising method to measure Ag 37, |39| ]. 

If parton kinematics can be approximately reconstructed, one can bin the 
events in the parton momentum fractions xi,X2 of the hard scattering. Assum- 
ing dominance of the Compton process, the asymmetry All for prompt-photon 
production can then be written as 



A 



LL 



Ag(xi) 



el [Ag(x2) + Ag(x2)] 
Ege2 [q{x2) + q{x2)] 



■ aLiigq ^ 79) + (1^2). (15) 



As a result of the quark charge-squared weighting, the second factor in Equa- 
tion |l^ coincides, to lowest order, with the spin asymmetry A^ measured in po- 
larized DIS, and the partonic asymmetry clll is calculable in perturbative QCD. 
Thus, from the measurement of All, one can directly extract Ag{x)/g{x). 

Both Phenix and Star intend to use this procedure for a direct leading-order 
determination of Ag, where one exploits the dominance of 2 — > 2 {ab — > 7c) parton 
scattering when reconstructing xi,X2- This is done either on average based on 
the detector acceptance for the photon, or event-by-event by observing photon- 
plus-jet events (Star). Estimates of the "background" from qq annihilation 
have been made [^]. Eventually, the aim will be a "global" QCD analysis of 
polarized prompt photon, and other RHIC and DIS, asymmetry data to determine 
the full set of polarized parton densities simultaneously, as is done routinely in 
the unpolarized case |41, |4^, 43 1. In this case, one can directly work from the 
spin asymmetries, and inclusion of, for instance, higher-order corrections is more 
readily possible. 

Figure |^ shows the level of accuracy Star can achieve in a direct measure- 
ment of Ag based on reconstructing parton kinematics in photon-plus-jet events. 
The solid lines show in each plot the input density employed for Ag{x), taken 



from |44|. The data points and the error bars show the reconstructed Ag{x) and 
its precision for standard luminosities in runs at ^/s = 200 GeV (open circles) 
and ^/s = 500 GeV (solid circles). 
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Figure 8: Sensitivity of Star measurements of Ag{x) in the channel pp 
7-Fjet -hX. 



High-pT^ photons can also be produced through a fragmentation process, in 
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which a parton, scattered or produced in a QCD reaction, fragments into a photon 
plus a number of hadrons. The need for introducing a fragmentation contribution 
is physicahy motivated by the fact that a QCD hard-scattering process may 
produce, again through a fragmentation process, a p meson that has the same 
quantum numbers as the photon and can thus convert into a photon, leading to 
the same signal. In addition, at higher orders, the perturbatively calculated direct 
component contains divergencies from configurations where a final-state quark 
becomes collinear to the photon. These singularities naturally introduce the need 
for nonperturbative fragmentation functions into which they can be absorbed. So 
far, the photon fragmentation functions are insufficiently known; information is 
emerging from the LEP experiments |45|. Note that all QCD partonic reactions 
contribute to the fragmentation component; thus, the benefit of having a priori 
only one partonic reaction (qq — > ^g) competing with the signal {qg — > 'jq) is lost, 
even though some of the subprocesses relevant to the fragmentation part at the 
same time result from a gluon initial state. Theoretical studies |46, H^, Hq, 4£] 
for photon production in unpolarized collisions, based on predictions 51] 
for the photon-fragmentation functions that are compatible with the sparse LEP 
data, indicate that the fragmentation component is in practice a small, albeit 
nonnegligible, effect. 

In the fixed-target regime, fragmentation photons are believed ||4^ to con- 
tribute at most 20% to the direct photon cross section. At collider energies, 
the fragmentation mechanism is estimated to produce about half of the observed 
photons; however, an "isolation" cut can be imposed on the photon signal in ex- 
periment. Isolation is an experimental necessity: In a hadronic environment, the 
study of photons in the final state is complicated by the abundance of vr'^s, which 
decay into pairs of 7s. If two photons are unambiguously detected in an event, 
their invariant mass can indicate whether they resulted from a vr'^ (or 7]) decay. 
However, either escape of one of the decay photons from the detector or merging 
of the two photons from vr" decay at high px fake a single photon event. The 
isolation cut reduces this background, since vr'^s are embedded in jets. If a given 
neighborhood of the photon is free of energetic hadron tracks, it is less likely that 
the observed photon came from vr'^ decay, and the event is kept; it is rejected 
otherwise. Traditionally, isolation is realized by drawing a cone of fixed aperture 
in ip-r] space around the photon [where (p is the photon's azimuthal angle and 
r] = — lntan(0/2) is its pseudorapidity, defined through its polar angle 6], and 
by restricting the hadronic transverse energy allowed in this cone to a certain 
small fraction of the photon transverse energy. In this way, the fragmentation 
contribution to single 7s, resulting from an essentially collinear process, will also 
be diminished |52]. It is not expected |p7| , |4^ ] that fragmentation will remain 
responsible for more than 15-20% of the photon signal after isolation. It has 
been suggested [53] that allowing proportionally less hadronic energy, the closer 
to the photon it is deposited, rather than permitting a fixed fraction in the full 
isolation cone, would improve isolation by reducing the fragmentation photons 
still further . 

Several early theoretical studies for isolated prompt-photon production at po- 
larized RHIC have been published (e.g. [36, |3^, |3^, |3^). The QCD corrections to 
the direct (i.e. nonfragmentation) component of polarized prompt-photon produc- 
tion were first calculated in References [E3, |55| and are now routinely included in 



theoretical studies (e.g. [pq, 57, 58, 59,pn[). In particular. References [pq, 5£ 



present Monte Carlo codes for the next-to- leading-order (NLO) corrections to the 
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direct part of the cross section, which allow the isolation constraints to be taken 
into account and also have the flexibility to predict photon-plus-jet observables, 
pp — > 7 + jet + X. We also emphasize that much effort has gone, and is still 
going, into event-generator studies |4C, 63] for prompt-photon physics at 

RHIC. 
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Figure 9: Asymmetry as a function of transverse momentum, for various po- 



larized parton densities, at different cms energies [30|. The expected statistical 
errors for the Phenix experiment are also shown. 



Figure ^ shows the asymmetry as obtained in an NLO theory calculation, 
as a function of the photon's transverse momentum pT- A rapidity cut \rj\ < 
0.35 has been applied, matching the acceptance of the Phenix experiment. In 
the left (right) part of the figure we plot the asymmetries obtained at \/s = 
200 GeV (500 GeV). The isolation of Reference ||5^ was used, with isolation 
cone opening Rq = 0.4 and e-y = l,n = 1 (see Reference |Q for details on 
the latter parameters). The solid, dashed, and dotted lines correspond to the 
NLO predictions obtained with GRSV-STD, GRSV-MAXg ||, and GS-C || 
polarized parton densities, respectively. These densities are all compatible with 
present data from polarized DIS and differ mainly in their gluon content: GRSV- 
MAXg has a very sizeable positive gluon distribution, whereas GS-C has a small, 
and oscillating, Ag. The gluon of GRSV-STD lies between the other two. The 



three gluon densities are shown in Figure 22 in the Appendix. The error bars 
represent the expected statistical accuracy for the measurement at Phenix, with 
A(j) = TT and for standard luminosities and beam polarizations. 

It is a striking feature of Figure |^ that different spin-dependent gluon densities 
do indeed lead to very different spin asymmetries for prompt-photon production. 
RHIC experiments will be able to measure Ag'. 

For fixed pT, higher-energy probes lower x in the parton distributions, and this 
leads to the smaller predicted asymmetries for ■y/s=500 GeV. If one considers the 
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same xt = 2pt/\^ value for the two energies in Figure ^ the parton densities 
are being probed at similar momentum fractions but rather different "resolution" 
scales, of the order of pt- It will be interesting to see whether measurements 
performed at different cms energies will yield information that is consistent, and 
compatible, with QCD evolution. 

Present comparisons between theory and experiment [and possibly between 
experiment and experiment [^]] regarding unpolarized direct 7 production are 
unsatisfactory |65]. Transverse smearing of the momenta of the initial partons 



participating in the hard scattering, substantially larger than that already in- 



troduced by the NLO calculation, has been considered |6£, 41, ^ to reconcile 
theory with data. This approach is partly based on measured values of dimuon, 
dijet, and diphoton pair transverse momenta in hadronic reactions [^] and 
has enjoyed some phenomenological success. More recently, the role of all-order- 
resummations of large logarithms in the partonic cross section, generated by (mul- 
tiple) soft-gluon emission, has been investigated in the context of prompt-photon 



production p8| , 69, 70, 71]. Threshold resummations have been shown 



to lead to improvements in the fixed-target regime, and a very recent new for- 



malism 1 71 1 that jointly incorporates threshold and /ct resummations has the 
potential of creating the substantial enhancements needed for bringing theory 
into agreement with data. It is likely that a better understanding of the prompt- 
photon process will have been achieved by the time RHIC performs the first 
measurements of pp ^ ^X. Also, the main problems reside in the fixed-target 
region; at colliders there is much less reason for concern. RHIC itself should also 
be able to provide new and complementary information in the unpolarized case — 
never before have prompt-photon data been taken in pp collisions at energies as 
high as =200-500 GeV. 



Finally, we note that it was also proposed |^, 72, 73 1 to determine lS.g through 
the reaction qg — > 7*^, which is again the Compton process, but now with a 
photon off-shell by the order of a few GeV and giving rise to a Drell-Yan lepton 
pair of comparable pT- The advantage is a cleaner theoretical description; for 
instance, no photon fragmentation component is present in this case. However, 
compared to prompt-photon production at a given pT, the event rate is reduced 
by 2-3 orders of magnitude due to the additional factor Oem/ (SvrQ^) in the Drell- 
Yan cross section, where Q is the dilepton mass. Higher statistics are available 
at lower pT, but at the price of reduced asymmetry and higher background from 
qq "f*g annihilation. 

3.2 Jet Production 

Toward the higher end of RHIC energies, jets could be the key to /S.g: at y/s = 
500 GeV, clearly structured jets will be copiously produced, and jet observables 



will show a strong sensitivity to lS.g thanks to the dominance |3£, 74] of the gg 
and qg initiated subprocesses (see Figure 06) in accessible kinematical ranges. Jet 
studies will be performed by Star. One can alternatively look for high-p^^ leading 
hadrons such as vr^ , vr^ , whose production proceeds through the same partonic 
subprocesses but involves an explicit fragmentation function in the theoretical 
description. This is planned for the Phenix experiment, where the limitation in 
angular coverage precludes jet studies. 

Knowledge of the NLO QCD corrections is expected to be particularly impor- 
tant for the case of jet production, since it is only at NLO that the QCD structure 
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of the jet starts to play a role in the theoretical description, providing for the 
first time the possibility of realistically matching the procedures used in experi- 
ment in order to group final-state particles into jets. The task of calculating the 



NLO QCD corrections to polarized jet production has been accomplished |75|. 
Furthermore, a Monte Carlo code that had been designed by Frixione ||76|, based 
on Reference [77| and the subtraction method in hadron-hadron unpolarized col- 



lisions, was extended to the polarized case in Reference We emphasize that 
in the unpolarized case, the comparison of NLO theory predictions with jet pro- 
duction data from the Tevatron is extremely successful (see e.g. |7^). 



Figure |10| shows the double-spin asymmetry for single-inclusive jet production 



at NLO as a function of the jet pT and for various polarized parton densities |44, 



79 1 with different Ag (see Figure 22 in the Appendix). A cut |r/| < 1 has 



been applied, and we have chosen the Ellis-Soper (ES) cluster jet algorithm |8C] 
with the resolution parameter D = 1. The renormalization and factorization 
scales have been chosen as //q ~ Pt (for further details, see Reference ||75|). The 
asymmetry shows a strong sensitivity to A^f. However, the asymmetry is rather 
small, regardless of the specific parton densities used. Fortunately, the expected 
statistical accuracy of such a jet measurement, calculated for standard luminosity 
and indicated in the figure, is very good. 
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Figure 10: Asymmetry versus jet transverse momentum [|75[ for various polarized 
parton density sets. The "data point" for pT=48 GeV/c indicates the statis- 
tical accuracy expected for the Star experiment for standard polarization and 
luminosity. Expected errors for lower pT are smaller than the points shown. 



The inclusion of NLO corrections in jet production, as shown in Figure 11 
leads to a clear reduction in scale dependence of the cross section. One thereby 
gains confidence that it is possible to calculate reliably the cross section and the 
spin asymmetry for a given Ag. This reduction in scale dependence after NLO 
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Figure 11: Scale dependence of the next-to- leading-order and Born pj- distri- 
butions for jet production |75|. (a) Polarized pp scattering and (6) unpolarized 
pp scattering at ^/s = 500 GeV. The range of the pseudorapidity is restricted to 
\v\ < 1- 



3.3 Heavy-Flavor Production 

The production of heavy quark pairs in hadronic collisions is dominated by gluon- 
gluon fusion, gg QQ (see Figure]^). For pp collisions, the competing channel 
qq — > QQ is particularly suppressed, since it requires an antiquark in the initial 
state. Thus, heavy quarks provide direct access to the gluons in the proton. 
Early predictions [^l| at the lowest order demonstrated that indeed this reaction 
could be used to measure Ag in polarized pp collisions. The importance of NLO 



corrections for a quantitative analysis was pointed out |82]. Presently, only the 



NLO QCD corrections to heavy-flavor production in polarized photon-photon 



and photon-proton [gj, 85 1 collisions are known; it is anticipated that the full set 
of NLO corrections relevant for polarized pp collisions will be available soon [^]. 
It should be mentioned that in the unpolarized case, theoretical NLO predictions 
for hadro- and photoproduction of heavy flavors often fail to provide a satisfactory 
description of the data (see [^7| for review) . 

Heavy-flavor production can be selected by the channels pp P^X, pp — > 
e^X, pp p^fi^X, pp — > e^e~X, and pp — > p^e^X. Like sign leptons are 
also possible from bottom, with one direct 6-decay to a lepton and one sequential 
decay through charm. Charm and bottom events will probe the gluon density at 
different momentum fractions and scales, and also enter the analysis with differ- 
ent, albeit calculable, weights. Experimentally it may be possible to determine 
the fraction of the charm production rate by, for example, looking at the channel 
pp p+D^X. 

The production of heavy quarkonia is another potentially attractive probe of 
the gluon density with a clear experimental signature. However, so far we do not 



24 



Bunce et al 



understand the production mechanism. Predictions for tp production based on the 
color-singlet model |88| fall short of experimental data taken at the Tevatron (see 
e.g. [^). This has stimulated the development of a more general approach that 
also gives rise to potentially important color-octet contributions [^]. Theoretical 
studies for the spin asymmetry in charmonium production in pp collisions have 
|91|, m, |9|, H, 11]. Reference |2| considers the color -singlet 
also examines color-octet contributions. Sensitivity 



been presented 
mechanism; Reference 




to the production mechanism as well as to Ag is found. Similarly, X2(3556) 
production at RHIC would have the potential to discriminate between the color- 
singlet and the color-octet mechanisms, as well as to pin down Ag [Q. Here, 
one would have to look at the angular distribution of the decay photon in X2 
J /ip + ^. The number of observed events for this reaction will unfortunately be 
low at RHIC. 
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4 QUARK AND ANTIQUARK HELICITY DISTRIBUTIONS 

Measurements in polarized DIS when combined with information from baryon 
octet /3-decays [Q], show that the total quark-plus-antiquark contribution to the 
proton's spin, summed over all flavors, is surprisingly small. In the standard in- 
terpretation of the /3-decays Q, this finding is equivalent to evidence for a large 
negative polarization of strange quarks in the proton, which makes it likely that 
also the SU{2) {u,d) sea is strongly negatively polarized. This view is corrobo- 
rated by the fact that in this analysis the spin carried, for example, by u quarks 
comes out much smaller than generally expected in quark models implying 
that a sizeable negative n-sea polarization partly compensates that of the valence 
u quarks. Alternative treatments of the information from /3-decays when 
combined with the DIS results, also directly yield large negative u and d polar- 
izations. Inclusive DIS (through 7* exchange) itself is sensitive to the combined 
contributions of quarks and antiquarks of each flavor but cannot provide informa- 
tion on the polarized quark and antiquark densities separately (see Appendix). 
Directly measuring the individual polarized antiquark distributions is therefore 
an exciting task and will also help to clarify the overall picture concerning DIS 
and the /3-decays. 

Further motivation for dedicated measurements of antiquark densities comes 
from unpolarized physics. Experiments in recent years have shown [9^ , |97| , a 
strong breaking of SU{2) symmetry in the antiquark sea, with the ratio d{x)/u{x) 
rising to 1.6 or higher. It is very attractive to learn whether the polarization of u 
and d is large and asymmetric as well. RHIC experiments will measure the d/u 
unpolarized ratio and the u and d polarizations separately. 

Semi-inclusive DIS measurements are one approach to achieving a separa- 
tion of quark and antiquark densities. This method combines information from 
proton and neutron (or deuteron) targets and uses correlations in the fragmen- 
tation process between the type of leading hadron and the flavor of its parton 
progenitor, expressed by fragmentation functions. The dependence on the de- 
tails of the fragmentation process limits the accuracy of this method. At RHIC 
the polarization of the u, u, d, and d quarks in the proton will be measured di- 
rectly and precisely using maximal parity violation for production of W bosons 



in ud and du W |Tog, |lOl], |0|, |To|. In addition, at RHIC, 



inclusive production of tt, K, and A will be used to measure quark and antiquark 
polarization through the fragmentation process. Another probe at RHIC will be 



Drell-Yan production of lepton pairs [104, M, 39, pa, |0|, no|, 103]. 



4--1 Weak Boson Production 

Within the standard model, W bosons are produced through pure V-A interac- 
tion. Thus, the helicity of the participating quark and antiquark are fixed in the 
reaction. In addition, the W couples to a weak charge that correlates directly 
to flavors, if we concentrate on one generation. Indeed the production of Ws in 
pp collisions is dominated by u,d,u, and d, with some contamination from s, c, s, 
and c, mostly through quark mixing. Therefore W production is an ideal tool to 
study the spin-flavor structure of the nucleon. 

The leading-order production of Ws, ud W~^ , is illustrated in Figure 
The longitudinally polarized proton at the top of each diagram collides with an 
unpolarized proton, producing a W~^ . At RHIC the polarized protons will be 
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Figure 12: Production of a in a, pp collision, at lowest order, (a) Au is probed ii 
the polarized proton. (6) Ad is probed. 



in bunches, alternately right- (+) and left- (— ) handed. The parity-violating 
asymmetry is the difference of left-handed and right-handed production of Ws, 
divided by the sum and normalized by the beam polarization: 



AY 



1 N^{W)-N+{W) 



P " N^{W) + N+{W) ■ ^^^^ 

As Figure ^ shows, we can construct this asymmetry from either polarized beam, 
and by summing over the helicity states of the other beam. The production of 
the left-handed weak bosons violates parity maximally. Therefore, if for example 
the production of the proceeded only through the diagram in Figure 
the parity-violating asymmetry would directly equal the longitudinal polarization 
asymmetry of the u quark in the proton: 



Similarly, for Figure \L 



AT 



u_{xi)d{x2) — Uj^{xi)d{x2) /^u{xi) 

uZ{xi)d{x2) + u^{xi)d{x2) u{xi) 

alone, 

d1:{xi)u{x2) - dXixi)u{x2) _ Ad{xi) 
d~t{xi)u{x2) — d'^{xi)u{x2) d{xi) 



(17) 



(18) 
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In general, the asymmetry is a superposition of the two cases: 



Ar 



Au{xi)d{x2) — Ad{xi)u{x2) 
u{xi)d{x2) + d{xi)u{x2) 



(19) 



To obtain the asymmetry for , one interchanges u and d. 

For the pp cohisions at RHIC with ^/s = 500 GeV, the quark wih be predom- 
inantly a valence quark. By identifying the rapidity of the W, yw, relative to 
the polarized proton, we can obtain direct measures of the quark and antiquark 
polarizations, separated by quark flavor: approaches Au/u in the limit of 

yw ^ 0, whereas for yw <C the asymmetry becomes —Ad/d. Higher-order 
corrections change the asymmetries only a little |102| , |103| ]. 

The kinematics of W production and Drell-Yan production of lepton pairs is 
the same. The momentum fraction carried by the quarks and antiquarks, xi and 
X2 (without yet assigning which is which), can be determined from yw, 



Xi 



Mw 



,yw 



X2 



Mw 



-yw 



(20) 



Note that this picture is valid for the predominant production of Ws at pT = 0. 
The experimental difficulty is that the W is observed through its leptonic decay 
W Iv, and only the charged lepton is observed. We therefore need to relate the 
lepton kinematics to yw, so that we can assign the probability that the polarized 
proton provided the quark or antiquark. Only then will we be able to translate 
the measured parity-violating asymmetry into a determination of the quark or 
antiquark polarization in the proton. 

The rapidity of the W is related to the lepton rapidity in the W rest frame 
{yD and in the lab frame {y\^^) by 



Vi"'' = Vi + yw, where = -In 



1 + cos6'* 
1 — cos6'* 



(21) 



Here 9* is the decay angle of the lepton in the W rest frame, and cos9* can be 
determined from the transverse momentum (pr) of the lepton with an irreducible 
uncertainty of the sign [107|, since 



lepton * 

Pt = Pt 



Mw 



sm9*. 



(22) 



In this reconstruction, the pt of the W is neglected. In reality, it has a pt, 
resulting for example from higher-order contributions such as gu — > W^d and 
ud W^g, or from primordial px of the initial partons. 

Usually W production is identified by requiring charged leptons with large pr 
and large missing transverse energy, due to the undetected neutrino. Since none 
of the detectors at RHIC is hermetic, measurement of missing pj- is not avail- 
able, which leads to some background. Possible sources of leptons with high px 
include charm, bottom, and vector boson production. Above px > 20 GeV/c, 
leptons from W decay dominate, with a smaller contribution from produc- 
tion. Both Phenix and Star can estimate the single-lepton background from 
measured Z^ production. The additional background from misidentified hadrons 
is expected to be small. 

Expected yields were estimated with Pythia yi08| and ResBos [ [L09| |. The 
cross section at RHIC for production is about 1.3 nb (0.4 nb). These 
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estimates vary by 5-10% according to the choice of the parton distribution set. 
For 800 pb-i and pr > 20 GeV/c, Phenix expects about 8000 W+s and 8000 
W~s in the muon arms (that the numbers are equal is due to the decay angle 
distribution and acceptance), as well as 15,000 and 2500 W~ electron decays 
in the central arms. Star, with its large acceptance for electrons, expects 72,000 
W~^s and 21,000 W~s. Using Equation ^to reconstruct x, Figure 13 shows the 
expected sensitivity for Af{x)/f{x), with / = u,d,u,d, for the Phenix muon 
data. 




Figure 13: Expected sensitivity for the flavor-decomp osed quark and antiquark polar- 
ization overlayed on the parton densities of Reference |101] (BS) and of Reference jiij 
[GS95LO(A)]. Darker points and error bars refer to the sensitivity from Al{W'^) mea- 
surements, and lighter ones correspond to Al{W~). 



RHIC will also significantly contribute to our knowledge about the unpolar- 
ized parton densities of the proton, since it will have the highest-energy pp col- 
lisions, pp production of Ws has a much stronger valence component in the 
determined |11C| ] u{x)/d{x) ratio. Isospin dependence in Drell-Yan production of 
muon pairs in pp,pd scattering [p^ j, violation of the Gottfried sum rule | 111 , 96 1, 
and recent semi-inclusive DIS measurements [98| have shown that the unpolar- 
ized sea is not SU{2) symmetric. At RHIC, the ratio of unpolarized and 
W~ cross sections will directly probe the d/u ratio, as shown in Figure 14. 
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Figure 14: The ratio Rw = {da{W~) / dy) / {da{W~^) / dy) for unpolarized pp collisions 
at RHIC. The shaded region indicates that unpolarized pp collisions are symmetric in 
yw- To illustrate the sensitivity of the measurement, we show an earlier set of parton 



densities (CTEQ4M |112|] ) and a set (MRS99 flij) that includes the latest information 
from Drell-Yan data |97| . Both curves include an asymmetric sea with d/u rising to 1.6 
for increasing antiquark momentum fraction Xq, but the latter also includes a drop-off in 
the ratio for higher Xq. 



4-.2 Drell-Yan Production of Lepton Pairs 

Drell-Yan production of lepton pairs has been a basis for information about sea 
quarks |^]. At lowest order, lepton pairs are created from quark-antiquark an- 
nihilation. With knowledge of the quark densities, Drell-Yan cross sections give 
the antiquark distributions versus x. The spin asymmetry An for Drell-Yan 
lepton pair production in collisions of longitudinally polarized proton beams is 
proportional to a sum of contributions over quark flavors, each a product of the 
polarized quark density times the antiquark distribution. The subprocess ana- 
lyzing power is maximally negative, aiL = —1- One therefore has, at lowest 
order. 

All = aLL x + Ag(xi)Ag(x2)} 

J2gej{qixi)q{x2) + q{xi)q{x2)} 

This asymmetry is parity-conserving if the process proceeds via a photon. Since 
the cross sections by flavor are weighted by the electric charge squared, the 
asymmetry is dominated by the uu combination and gives information on the 
u polarization, with the u quark polarization as input. NLO corrections to the 
asymmetry have been calculated [|106|, |102|] to be small for low px of the vir- 
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tual photon. For higher pT, Dreh-Yan production is sensitive to Ag(x) through 
Qg ^ 1*1 1 37, 72, 73], as discussed in Section ^. 

However, lepton pair production in high-energy pp cohisions is dominated by 
coincidental semileptonic decays of heavy-quark pairs, e.g. 6 — > cZ~ in the low- 
mass region. The feasibility of the measurements will therefore depend on the 
ability to separate or estimate this background. Estimates of the yields in the 
Phenix muon arms obtained with Pythia for pp collisions at y/s =200 GeV show 
that lepton pairs with invariant mass M >6 GeV/c^ are dominated by Drell-Yan 
production. One expects ~40,000 pairs for a nominal integrated luminosity of 
320 pb^i. 
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5 TRANSVERSE AND FINAL-STATE SPIN EFFECTS 

Exciting physics prospects also arise for transverse polarization of the RHIC pro- 
ton beams. One is the possibility of a first measurement of the quark transversity 
densities introduced in Table [l|. The transversity distributions, measuring dif- 
ferences of probabilities for finding quarks with transverse spin aligned and anti- 
aligned with the transverse nucleon spin, are as fundamental as the longitudinally 
polarized densities for quarks and gluons, Aq, Ag; they have evaded measurement 
so far because they decouple from inclusive DIS. Comparisons of the polarized 
quark distributions 5q and Aq are particularly interesting; in the nonrelativistic 
limit, where boosts and rotations commute, one has 6q{x,Q^) = Aq{x,Q'^). De- 
viations from this provide a measure of the relativistic nature of quarks inside 
the nucleon. 

Studies of single-transverse spin asymmetries, defined similarly to Equation 
will be a further interesting application. They arise as "higher-twist" effects (that 
is, they are suppressed by inverse powers of the hard scale) and probe quark-gluon 
correlations in the nucleon. They have an exciting history in experiments that 
were carried out at energies much lower than RHIC's, where large polarizations 



and single-spin asymmetries have been seen |114]. Yet another field of spin physics 
to be thoroughly examined by the RHIC experiments will be the transfer of 
longitudinal or transverse polarization from the initial into the final state, which 
then leaves traces in the polarization of hadrons produced in the fragmentation 
process. 

5. 1 The Quark Transversity Distributions 

The transversity densities 6q and 6q are virtually inaccessible in inclusive DIS 



17|. We can see this as follows [115]. In a simple parton model, and working in 
a helicity basis, we can view the quark densities as imaginary parts of polarized 
quark-hadron forward scattering in the ti-channel, denoted by A{H, h; H' , h') (see 
Figure One then has q = A{++; ++) + A{+-] +-), Aq = A{++; ++) - 

A{-\ — ;H — ), but 5q = ^(-|-+; ). Thus, for transversity to contribute, the 

quark has to undergo a helicity fiip in the hard scattering, which is not allowed 
(for massless quarks) at the DIS quark-photon vertex due to helicity conservation. 

Note the striking feature that the helicity labels of the final state in ^(+-|-; ) 

differ from those of the initial state. In other words, the complex conjugate 

amplitude contained in ^(+-|-; ) refers to a different physical state than the 

initial. This "off-diagonal" nature in terms of helicity is usually referred to as 



chiral-odd |11] and can indeed in practice only be achieved by having transverse 
polarization, which can be written as a superposition of helicity states. 

Another important consequence is that, unlike the situation for unpolarized and 
longitudinally polarized densities, there is no transversity gluon distribution [pl| . 



16 , 17]. This is due to angular momentum conservation; a gluonic helicity-flip 
amplitude would require the hadron to absorb two units of helicity, which a spin- 
1/2 target cannot do. 

The joint description of the quark distributions in terms of the A{H, h; H' , h') 
implies that transversity is not entirely unrelated to the q,Aq. Indeed, rewrit- 
ing ]|115| ] A{H,h;H',h') = J2x ^*H'h'(-^)^Hh{^)j where X is an arbitrary final 
state, one finds from the condition J2x l^++(^) ^ ^ — inequal- 
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ity ra 



q{x) + Aq{x) > 2\6q{x) 



(24) 



Figure 16 displays the region allowed by Equation which is indeed smaller 
than the one resulting from the trivial condition |(^(7(x)| < q{x). Equation 24 



holds for all quark flavors and separately for their corresponding antiquarks. As 
was demonstrated in References [118, |119|, 120, 121|, the inequality is preserved 



under QCD evolution; that is, if it is assumed to be satisfied at one resolution 



scale, it will hold at all larger scales. This remains true [119, 12C, 121| even at 
two-loop order [|119| , |122|| in evolution. 

The helicity flip required for transversity to contribute to hard scattering can 
occur if there are two soft hadronic vertices in the process. In this case, transverse 
spin can be carried from one hadron to the other along a quark line. One possi- 
bility is to have two transversely polarized hadrons in the initial state, as realized 
at RHIC. An alternative is to have one transversely polarized initial hadron and 
a final-state fragmentation process that is sensitive to transverse polarization. 
Here, the other initial particle could be a lepton, as in DIS, or another proton, 
as at RHIC. 

For the first possibility, a promising candidate process for a measurement of 
the Sq, 5q is Drell-Yan dimuon production which, to lowest order in QCD, pro- 
ceeds via qq 7* annihilation. A systematic study of this process was in fact 
also the place where the transversity densities made their first appearance in 



theory [15|. On the downside of this reaction is that the transversity antiquark 
density in the nucleon is presumably rather small; there is no splitting term 
g ^ qq in the evolution equations for transversity |l^, so a vital source for 
the generation of antiquarks is missing (only higher orders in evolution produce 
antiquarks carrying transversity |119| , |121| ). Also, in Drell-Yan, the event rate 
is generally low. However, when compared to other conceivable reactions in pp 
collisions that serve to determine parton densities, the Drell-Yan process has the 
advantage that to lowest order there is no partonic subprocess that involves a 
gluon in the initial state. If a reaction does have a gluon-initiated subprocess, its 



transverse double-spin asymmetry is expected to be suppressed |123| , |20[|. This 
is because gluons usually strongly contribute to the unpolarized cross sections in 
the denominator of the asymmetry, whereas they are absent for transversity, as 
discussed above. In addition, for many reactions other than Drell-Yan, one finds 



a particular "selection-rule" [ 123 , 2C | suppression of the contributing transverse 
subprocess asymmetries. 

Several phenomenological studies of Drell-Yan dimuon production at RHIC 





Figure 15: Left: quark densities as related to polarized quark-hadron forward scattering 
in the u-channel. Labels refer to helicities. Right: decoupling of transversity from deep- 
inelastic scattering |116|. Quark chirality is not changed by coupling to a photon or a 
gluon. 
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have been presented [|T2|, |T2|, jTO^, |T2|, |T2|, ^^]. Model estimates of the 



transversity densities have been obtained in these studies by either assuming 



25q{x,QQ) = Qo) + ^Qi^, Qo) (see Equation or by employing ||12^, |12|] 
6q{x,QQ) ~ Aq^XjQq), at some initial (low) resolution scale Qo- Note that the 
latter ansatz violates inequality 24 if Aq{x,QQ) < —^q{x,QQ). The transverse 
double-spin asymmetry for Drell-Yan dimuon production is (to lowest order) 



A' 



XT — CLXT' 



ge^g6q{xi,M^)6q{x2,M' 



+ (1^2) 



Ege2(;(xi,M2)^(x2,M2) + (1^2) 



(25) 



Here utt is the partonic transverse-spin asymmetry, calculable in perturbative 
QCD, and M is the dilepton mass. NLO corrections to Drell-Yan dimuon pro- 



duction with transversely polarized beams have been calculated [124, 125 , 121 , 
129 , |130(| and are routinely used in numerical studies. 

The Phenix endcaps will be able to identify muons with rapidity 1.2 < |y^j± | < 



2.4. Figure 17 shows predictions 121 | for Att- In order to model the transver- 
sity densities, saturation of inequality 24 at a low scale Q « 0.6 GeV has been 
assumed, making use of the information on the Aq, Aq in that inequality coming 
from polarized DIS. The statistical errors expected for Phenix are also shown. 
One observes that the asymmetry is generally small but could be visible experi- 
mentally if the transversity densities are not much smaller than those used here. 
Larger estimates for Att have been obtained |101], based on more optimistic 
assumptions concerning the size of the 5q, 5q. Careful studies of the background 
to lepton pair production resulting from coincidental semileptonic heavy-flavor 
decays (see Section Q) will be important. 

The other possibility involves one transversely polarized initial hadron and 
a final-state fragmentation process that is sensitive to transverse polarization. 
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Figure 17: Next-to-leading-order transverse-spin asymmetry for Drell-Yan dimuon pro- 
duction at = 200 GeV p^. 



Promising approaches have emerged from considering the production of high-p^^ 



dimeson systems | 131 , 12, 132], or from taking into account "intrinsic" transverse 



momentum degrees of freedom in a fragmentation process producing a single 



high-pT pion |133]. Both dimesons and pions are very abundantly produced in 



high-energy pp collisions. It has been shown [131] that the azimuthal distribu- 
tion of low-mass pairs of pions about the final-state jet axis can be used as a 
measure of the transverse polarization of the quark initiating the jet. The same 



is true [133] for the "intrinsic" transverse momentum distribution of a produced 
pion relative to its quark progenitor. In this way, one effectively obtains an asym- 
metry that is sensitive to products of the transversity density for the initial-state 
quark and a transverse-polarization-dependent fragmentation function for the fi- 



nal state. For instance, for the mechanism proposed for DIS in Reference [133], 
the fragmentation function would be 

Hi{z, k±) oc D^/g^ {z, k±) - D^/gi {z, k±) , (26) 

where k± is the "intrinsic" transverse momentum in the fragmentation process. 
Notice that one polarized proton in the initial state is sufficient for this kind of 
measurement. Time-reversal invariance, however, precludes a nonzero effect un- 
less phases are generated by final-state interactions in the fragmentation process 
that do not average to zero upon summation over unobserved hadrons. It is a pri- 
ori unclear whether such a net phase will exist. This led to investigation [132[ of 



the interference between s and p waves of two-pion systems with invariant mass 
around the p. Such an interference effect yields sensitivity to the polarization 
of the quark progenitor through the quantity k^+ x k^~ ■ st, where the ks are 
the pion momenta and st is the transverse nucleon spin; one effectively uses the 
angular momentum of the two-pion system as a probe of the quark's polarization. 
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Staying in the mass region around the p ensures ]132| ] that the final-state inter- 
action phase does not average to zero. The s-p wave interference in the q — > vrvr 
formation is described by a new set of fragmentation functions, the interference 
fragmentation functions |132|. Just as the function in Equation 26, the latter 
are presently entirely unknown; the price to be paid for obtaining sensitivity to 
transversity in all of the ways suggested in Reference [131, 132, 



331 is thus the 



introduction of another unknown component. However, one may hope that the 
involved fragmentation functions can be determined independently in e~^e~ anni- 
hilation. Studies of the experimental situation at RHIC concerning the proposal 



of 1 132 1 are under way |13 



5.2 Transverse Single-Spin Asymmetries 



Surprisingly large single-transverse spin asymmetries, for instance in fixed-target 
p^p — > nX at pion transverse momenta of a few GeV, have been observed ex- 
perimentally [114] over many years. RHIC will further investigate the origin of 
such asymmetries. Within the "normal" framework of perturbative QCD and the 
factorization theorem at twist-2 for collinear massless parton configurations, no 
single-transverse spin asymmetry is obtained — nonzero effects occur only when 
one keeps quark mass terms (as is required to generate helicity flips) and when 
one takes into account at the same time higher-order loop diagrams that pro- 



duce relative phases [|8|. Such effects are therefore of the order of agmq/y/s 
and cannot explain data such as that in Reference [114]. It is believed that 
nontrivial higher-twist effects are responsible for the observed single-spin asym- 

Reference 



metrics [135, 136, |13^ [. Reference |136| , |13^ ] showed how single transverse-spin 
asymmetries can be evaluated consistently in terms of a generalized factorization 
theorem in perturbative QCD, wherein they arise, for example, as convolutions 
of hard-scattering functions with an ordinary twist-2 parton density from the un- 
polarized hadron and a twist-3 quark-gluon correlation function representing the 
polarized hadron. Another contribution involves the transversity distribution and 



another (chiral-odd) spin-independent twist-3 function of the proton [138 
A simple model was constructed [113 



13C]. 



138 1 that assumes only correlations of va- 



lence quarks and soft gluons. It can describe the present data and makes various 
definite predictions, to be tested at RHIC, where one certainly expects to be in 
the perturbative domain. In particular, at RHIC, one should see the fall-off with 
Pt of the single-transverse spin asymmetries in single-inclusive pion production, 
associated with their twist-3 nature (see Figure 18). 

A related dynamical origin for transverse single-spin asymmetries was pro- 
posed 133, 14C] to reside in the dependences of parton distribution and 
fragmentation functions on intrinsic parton transverse momentum k^- In fact, 
the proposal of ]|133[] for measuring transversity in the proton, which we dis- 
cussed in the previous subsection, proceeds for pp scattering exactly through a 
single-transverse spin asymmetry, making use of the /cT-dependent fragmentation 
function in Equation Suppression of the asymmetry should also arise here, 
through a factor {kT)/pT- It has also been considered that single-spin asymme- 
tries might be generated by dependences of the parton distribution functions 
in the initial state ]p!9| , 140]. Here, one could have 



flT{x,k_L) 

hi{x, k±) 



/g/pT {X,k±)- fg/pl {X,k±) , 
fgypix, k±) - fql/pix, k±) 



(27) 
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Figure 18: (a) Experimental data |137| and theoretical calculations |138| for transverse 
single-spin asymmetries for 7r+ and tt^ production in pp collisions at -^5=20 GeV as 
functions oixp- Predictions for RHIC for pt=4 GeV/c are superimposed. The transverse 
momentum dependence for RHIC at 2;^? = 0.4 is shown in (b). 



as the driving forces. There is a qualitative difference between the functions in 
Equations ^ and 26: In order to be able to produce an effect, the latter requires 
final-state interactions (which are certainly present), to make the overall process 
time-reversal-symmetry-conserving (see the previous subsection). In contrast, 
the distributions in Equation ^ rely on the presence of nontrivial (factorization- 
breaking) initial-state interactions between the incoming hadrons [|141[| , or on 
finite-size effects for the hadrons | 142 |; they vanish if the initial hadrons are 
described by plane waves. This makes the "Collins function" (Equation 26) 
perhaps a more likely source for single-spin asymmetries. The reservations con- 
cerning Equation 27 notwithstanding, when a factorized hard-scattering model is 
evoked, each mechanism described by Equations 26 and ^ can by itself account 
the present p^p 



for 1 141, |140| , |143| ] the present p'p — > ttX data. Also, all could be at work simul- 
taneously and compete with one another. Single-spin Drell-Yan measurements 
at RHIC should be a good testing ground [14C] for the existence of effects related 
to Equation 27, since for Drell-Yan the Collins function (Equation 26) cannot 
contribute. 



5.3 Spin- Dependent Fragmentation Functions 

Even in the context of a parity-conserving theory like QCD, an asymmetry can 
arise for only one longitudinally polarized particle in the initial state, if the longi- 
tudinal polarization of a particle in the final state is observed. The measurement 
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of the polarization of an outgoing highly energetic particle certainly provides 
a challenge to experiment. A baryons are particularly suited for such studies, 
thanks to the self-analyzing properties of their dominant weak decay, A pvr" . 
Recent results on A production reported from LEP [[144j| have demonstrated the 
feasibility of successfully reconstructing the A polarization. 

Spin-transfer asymmetries give information on yet unexplored spin effects in the 
fragmentation process. For our A example, the longitudinal transfer asymmetry 
will be sensitive to the functions 

ADHz) ^ Df^^^\z) - Df^':;^\z) (28) 

describing the fragmentation of a longitudinally polarized parton i = q,q,g into 
a longitudinally polarized A, where D^^^(z) {Df^^^\z)) is the probability of 
finding a A with positive (negative) helicity in a parton i with positive helicity, 
carrying a fraction z of the parent parton's momentum (see Section |2|). As was 
shown in Reference |145, 146], the LEP measurements [ p.44]] have provided initial 



information on some combinations of the AD^ but leave room for very different 
pictures of the spin-dependence in A fragmentation. Measurements of the polar- 



ization of As produced in pp collisions at RHIC should vastly improve |147] , |l4g| 1 
our knowledge of the AD^. Figure 19 illustrates this by showing the longitudinal 



spin transfer asymmetry at RHIC, defined in analogy with Equation ^ as 

where the lower helicity index refers to the polarized proton and the upper to 
the produced A. Various models for the AD^, all compatible with the LEP 



data, have been used in Figure 19. It will be interesting to see which scenario is 



favored by the RHIC measurements. A cut of xt > 0.05 has been applied in the 



figure. As are very copiously produced at RHIC [147], resulting in small expected 
statistical errors. 



Similarly optimistic conclusions have been reached [15C] for the case of trans- 



verse polarization of one initial beam and the A, in which case RHIC experiments 
would yield information on the product of the proton's transversity densities and 
the transversity fragmentation functions of the A, which are both so far unknown. 
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Figure 19: Th e longitudinal spin transfer asymmetry in A production at RHIC {^/s 
500 GeV) 



147 1, as a function of rapidity of the A for various sets of spin-dependent 



fragmentation functions proposed in Reference [145|. For scenario 1, only strange quarks 
transmit polarization to the A. In scenario 2, there is also a (negative) contribution from 
up and down quarks [ |149| . In scenario 3, all quarks equally produce polarized As. The 
expected errors for Star with standard luminosity and polarization are comparable to 
the "data" shown for 77 = ±2, and smaller for the other points. 
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6 PHYSICS BEYOND THE STANDARD MODEL 



So far we have discussed probing the proton spin structure at RHIC, and both 
using and testing perturbative QCD in the spin sector. Spin is also an excehent 
tool to go beyond the standard model and to uncover important new physics, 
if it exists. Many extensions of the standard model have been proposed. Our 
purpose in this section is to illustrate this new potentiality by means of a specific 
example. 

Let us consider one-jet inclusive production. As discussed in Section 3, the 
cross section is dominated by the pure QCD gg, gq, and qq scatterings, but the 
existence of the electroweak interaction, via the effects of the and Z gauge 
bosons, adds a small contribution. Consequently, the parity-violating helicity 



asymmetry Al, defined as [|151 | 



A, 



da 



pp-^jetX 



cLEt 



dEr 



+ 



dEr 



-1 



(30) 



is expected to be nonzero from the QCD-electroweak interference (as shown in 

Mw,z/'^ is seen, which is the 
The cross sections are 



Figure 20). Additionally, a small peak near E'j' 
main signature of the purely electroweak contribution, 
for one longitudinally polarized beam, colliding with an unpolarized beam. The 
existence of new parity-violating interactions could lead to large modifications of 
this standard-model prediction [151]. 

First let us recall that the sensitivity to the presence of some new quark- 
quark contact interactions has been analyzed in Reference [152]. Such a contact 
interaction could represent the effects of quark compositeness, under the form 



r - 

t-qqqq — £ 



(31) 



where ^' is a quark doublet, A is a compositeness scale, and e = ±1. If parity 
is maximally violated, r/ = ±1. Figure 2C shows how the standard-model pre- 
diction will be affected by such a new interaction, assuming A = 2 TeV, which 
is close to the present limit obtained for example by the D0 experiment at the 
Tevatron [153[. The statistical errors shown are for standard RHIC luminosity 
of 800 pb^^, and for jets with rapidity \y\ <0.5, and include measuring A^ us- 
ing each beam, summing over the spin states of the other beam. Due to the 
parity-violating signal's sensitivity to new physics, RHIC is surprisingly sensitive 
to quark substructure at the 2- TeV scale and is competitive with the Tevatron, 
despite the different energy ranges of these machines. Indeed, a parity-violating 
signal beyond the standard model at RHIC would definitively indicate the pres- 
ence of new physics 



151[[ . 



RHIC-Spin would also be sensitive to possible new neutral gauge bosons [154]. 
A class of models, called leptophobic Z' , is poorly constrained up to now. Such 
models appear naturally in several string-derived models [ p.55| (nonsupersym- 
metric models may be also constructed [156[). In addition, in the framework 
of supersymmetric models with an additional Abelian gauge factor f/(l)', it has 
been shown |157| [ that the Z' boson could appear with a relatively low mass 
{Mz < Mz' <1 TeV) and a mixing angle with the standard Z close to zero. 
The effects of different representative models are shown in Figure pOl (see Refer- 



ence |154[ for details). RHIC covers some regions in the parameter space of the 



40 



Bunce et al 



different models that are unconstrained by present and forthcoming experiments 
(e.g. Tevatron Run II), and RHIC would also uniquely obtain information on the 
chiral structure of the new interaction. 
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Figure 20: Al for one-jet inclusive production in pp collisions, versus transverse en- 
ergy, for y/s — 500 GeV. The solid curve with error bars represents the standard-model 
expectation. The error bars show the sensitivity at RHIC for 800 pb~^, for the Star 
detector. The other solid curves, labeled by the product of erj, correspond to the contact 
interaction at A = 2 TeV |15S]. The dashed and dotted curves correspond to different 
leptophobic Z' models |f54|. The calculations are at leading order. 



Other possible signatures of new physics at RHIC have been investigated. Par- 
ticularly interesting quantities |10Cl| , |158| , |159(| are transverse (single or double) 
spin asymmetries for production, since these are expected to be extremely 
small in the standard model |IO0| , |I60(| . For instance, the case of the corre- 
sponding standard-model double spin asymmetry A^j, was examined in detail 



recently |I60|. Non- vanishing contributions could arise here for example in the 
form of higher-twist terms, which would be suppressed as powers of M^/M^, 
where M is a hadronic mass scale and Mw the W mass. Other possible contri- 



butions were demonstrated in |16C] to be negligible as well. By similar arguments. 



also the corresponding single-transverse spin asymmetry for production, Afj, 



is expected to be extremely small in the standard model |159| ] . New physics ef- 
fects, on the contrary, might generate asymmetries at leading twist, for example 
through non-(y — ^) (axial)vector couplings of quarks to the W, or through tensor 
or (pseudo)scalar couplings, all of which would also have to violate CP in order 
to generate a single-spin asymmetry A^. In particular the latter asymmetry has 



been examined with respect to sensitivity to new physics effects at RHIC p9| . 
For a case study, the minimal supersymmetric extension of the standard model, 
with R-parity violation, was employed, which contains scalar quark-H^ interac- 
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tions and complex phases, resulting in CP- violating effects. The results of |159| ] 
show that in this particular extension of the standard model, is likely to be 
very small as well, below the detection limit of RHIC. Nevertheless, this does not 
exclude that other non-standard mechanisms produce larger effects, and and 
A^rp will be measured at RHIC with transversely polarized beams in the context 
of the physics discussed in the previous section. A non-zero result would be a 
direct indication of new physics. 
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7 SMALL-ANGLE pp ELASTIC SCATTERING 

In previous sections, we have discussed the physics of hard scattering at RHIC 
with polarized protons, which can be understood as colhsions of polarized quarks 
and gluons. The scattering is so energetic that we can use perturbative QCD 
to describe the interactions of the quarks and gluons, and, thus, probe the 
spin structure of the proton at very small distances. For example, scattering 
at Q2=(80 GeV)2 probes wave lengths of 0.003 fermi. Small-angle scattering, 
from total cross section to t = — 1 (GeV/c)^, probes the static proton properties 
and constituent quark structure of the proton, covering distances from 4 fermi 
[—t = 0.003 (GeV/c)^ in the Coulomb nuclear interference (CNI) region] to a 
distance of ~0.2 fermi. Unpolarized scattering shows striking behavior in this 
region, from the surprise that total cross sections rise at high energy, to observed 
dips in elastic cross sections around —t = l (GeV/c)^. The pp2pp experiment at 
RHIC will explore this region for spin-dependent cross sections, for ^/s=20- 
500 GeV, for the first time. 

Historically, new spin-dependent data have often shown new structure under- 
lying spin-independent cross sections, indicating the presence of unexpected dy- 
namics in the interaction. Several examples have been discussed in previous 
sections. Previous work with spin stops at y^=20 GeV, where tertiary polarized 
p and p beams were collected from the parity-violating decays of A and A hy- 
perons and steered onto unpolarized hydrogen and polarized pentanol {C5H12O) 
targets [161|. RHIC will provide much higher intensity, a large extension of the 



energy range, and pure targets for 2-spin measurements. 

In the energy regime ^/s >20 GeV, total cross sections have been observed 
to rise with energy for pp, pp, ir^p, and K^p. The pp total cross section rises 
through the Tevatron maximum energy of 2 TeV, and the pp total cross section 
has been observed to rise through its highest energy measurement at the ISR, 



y^=62 GeV [162]. The pp2pp experiment will measure spin-dependent total 
cross sections, (J||, a^i, and ai = a+ — a- [where the arrows represent transverse 
spin measurements, and (+) and (— ) represent helicities] through the range of 
rising cross sections available at RHIC. The unpolarized pp total cross section 
measurements will also be extended to -^5=500 GeV. 

For pp, the rise of the total cross section has been successfully described in the 
impact picture approach on the basis of the high-energy behavior of a relativistic 
quantum field theory |163| ] . This is based on the fact that the effective interaction 
strength increases with energy in the form s^~^'^/{ln sY , a simple expression in 
two key parameters c and c', where s is expressed in GeV^. A fit of the data 
then leads to the values of the two free parameters c = 0.167, c' = 0.748 1164] 



If this picture is correct (the field theoretical argument is based on connecting 
QED and QCD theories, but successfully predicted that the pp total cross section 
would continue to rise, following these parameters), there should be no difference 
in the rise of pp and pp total cross sections. An extension of this approach allows 
a description of the elastic cross section [165[, which will also be measured at 
RHIC. 

The single-spin asymmetry for pp elastic scattering, A^r, is expected to be 
small but significant in the CNI region, -t = 0.001-0.01 (GeV/c)^ As 
discussed previously, pp elastic scattering in the CNI region will be the basis 
of the RHIC polarimetry. CNI scattering is expected to produce an asymmetry 
from scattering an unpolarized proton (polarization averaged to zero) in one beam 
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from the magnetic moment of a polarized proton from the other beam, with a 
maximum of A]\i = 0.04 at —t = 0.003 (GeV/c)^. However, a hadronic spin- 
flip term can also contribute to the maximum, and this term is sensitive to the 



static constituent quark structure of the proton. The authors of Reference |16(:] 
remark that the helicity flip probes the shortest interquark distance in the proton, 
and that the helicity nonflip is sensitive to the largest quark separation in the 
proton due to color screening. The helicity-flip term, if present, can indicate an 
isoscalar anomalous magnetic moment of the nucleons [167], an anomalous color- 
magnetic moment causing helicity nonconservation at the constituent quark-gluon 
vertex |l68f| , and/or a compact quark pair in the proton |24, 16S ]. 



The only measurement of An in the CNI region at higher energy is by E704 
at Fermilab [161] at a lab momentum pl = 200 GeV/c; the results are shown 
in Figure m. The errors are too large to allow an unambiguous theoretical 
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Figure 21: Transverse single-spin asymmetry for proton proton elastic scattering. The 



data points are from Fermilab E704 [161|. The solid curve is the best fit with the spin- 
flip hadronic amplitude constrained to be in phase with non-flip hadronic amplitude; the 
dotted curve is the best fit without this constraint. 

interpretation. There are two fits to the E704 data shown with a nonzero 
hadronic spin flip term |2^]. As emphasized in References ]|2^ and 126], a large 



value of the hadronic helicity-flip amplitude generates a very large change in 
the maximum in Ajsf, which can be of the order of 30% or more. The pp2pp 
experiment will measure A]\f to ±0.001 in the CNI peak. This level of preci- 
sion is required for absolute polarimetry, giving an expected precision on 
of A^Ar/AAr=±0.001/0.04=±0.025. This experiment wih cover from 0.0005 < 
—t < 1.5 (GeV/c)^ (with additional detectors for the larger —t region). Thus, 
the location of the maximum in and its maximum value and shape will be 
determined. 

Small-angle scattering at high energy is presently understood in the Regge 
picture as being dominated by Pomeron exchange |17C| ]. The Pomeron, which 
has the vacuum quantum numbers with charge-conjugation C = +1, can be 
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interpreted as a two-gluon exchange. There is room in the data for a smah three- 
gluon exchange contribution with C = —1, the Odderon |171]. It has been shown 
recently |172| that the behavior of the two-spin transverse asymmetry ^atat in pp 
elastic scattering in the CNI region depends strongly on the Odderon contribution 
and that the pp2pp experiment is quite sensitive to its presence. 

In addition to the measurements discussed above, the pp2pp experiment will 
measure larger angles, to —t = 1.5 (GeV/c)^, which includes the region of dip 
structure in the unpolarized cross section, measuring A^- and the two-spin asym- 
metries Ann, Ass, and All \l§§,\- These first, precise, determinations of spin 
dependence for small-angle pp elastic scattering in the energy range = 20- 
500 GeV probe the spin structure of the proton in the nonperturbative QCD 
region, from the static properties of the proton to its constituent quark struc- 
ture. 

At higher energy, such as at the LHC, the CNI region becomes inaccessible. The 
minimum —t reachable with colliding beams depends on scattering the protons 
out of the beams. For fixed —t, the scattering angle falls as l/pbeam; whereas the 
beam size falls more slowly as 1/ ^/pbeam- Roughly, this limits an experiment at 
the LHC to -t >0.01 (GeV/c)^. 



8 CONCLUDING REMARKS 



RHIC will be the first machine to look at the proton spin structure by colliding 
polarized proton beams rather than scattering polarized leptons off polarized 
targets. Thus, one can test fundamental interactions in an entirely different 
environment and at much higher energies, as in the unpolarized case. (Here, 
too, information on the nucleon structure from DIS has been complemented by 
information from hadron colliders.) For hadron colliders, including RHIC-Spin, 
due to the high energy and luminosity that give access to hard parton scattering, 
perturbative QCD probes in one proton are used to study the nonperturbative 
structure of the "target" proton. 

What can we expect from RHIC-Spin? If, for example, a large gluon polariza- 
tion is observed, such a signal would imply a previously unknown fundamental 
role of the gluons in the proton spin. Surprise and new insights are very likely. 

This field is very new both theoretically and experimentally. Previous experi- 
mental spin work with hadron probes was at much lower energy and luminosity, 
and used impure polarized targets. Much of the discussion presented here is from 
very recent work. Thus, this article should not be seen as a review but rather as 
an invitation. 
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Appendix: Information from Polarized Deep-Inelastic Scattering 

In this Appendix we briefly discuss the information from DIS on Aq, Aq, Ag. If 
we neglect contributions resulting from or exchange, DIS is sensitive only 
to the sums of quarks and antiquarks for each flavor. Therefore, we define 

AQ{x,Q'^) = Aq{x,Q'^) + Aq{x,Q^). (32) 

To lowest order, we can then write the structure functions , 5" appearing in 
DIS off polarized proton and neutron targets as 

2gP{x,Q^) = ^AU{x,Q^) + 1\aV{x,Q^) + AS{x,Q^] 
9 9 L 

2g'^{x,Q^) = ^AV{x,Q^) + ^[aU{x,Q^) + AS{x,Q^)\ , (33) 

where all parton densities refer to the proton. We can compactly rewrite this as 

<7f "(x, Q2) = ±1a^3(x, + ^AAsix, + iAS(x, Q2) , (34) 

where the upper sign refers to the proton, and where we have introduced the 
flavor-non-singlet combinations A^s = AU - AV, A As = AU + AV - 2 AS, 
and the singlet AE = AU + AT> + AS. Had we data at only one Q^, the two 
structure functions Qi'"' could not provide enough information to determine the 
full set A^3, A^8, AS at this Q^. When information at different is available, 
one can combine the data with knowledge about QCD evolution. In particular, 
each non-singlet quantity evolves separately from all other quantities, whereas AS 
mixes with the polarized gluon density Ag{x, Q"^) in terms of a matrix evolution 



equation [10, |173|]. Thanks to this property under evolution, g\' {x,Q ) give in 



principle access to all four quantities, A^3,A^8,AS, and Ag [174, |175|] . We 



note that, when performing fits to data in practice, one usually also includes 
constraints on the "first moments" (Bjorken-x integrals) of A^3^8 derived from 
the /3-decays of the baryon octet, the constraint on A^3 being essentially the 
Bjorken sum rule |176| ]. In this way, one is also able to better determine the first 
moment of AS, which corresponds to the fraction of the proton spin carried by 
quarks and antiquarks. 

Information on A^3, A^g; and Ag is equivalent in a "three- flavor world" 
to information on AZY,AP,A5, and Ag — this is what DIS data can provide in 
principle. We emphasize again that inclusive DIS cannot give information on the 
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quark and antiquark densities separately; it always determines only the AQ. To 
distinguish quarks from antiquarks, let alone to achieve a full flavor separation 
of the polarized sea, one needs to defer to other processes (see Section |^. 

We do not address in detail the question of how well the present data, within 
their accuracies, do indeed constrain the quantities A^3, A^g; AS, and Ag. For 
this we refer the reader to the growing number of phenomenological analyses 
of the polarized DIS data pi, |4 



174, 79, 17 



|175| ]. However, to give a very 
rough picture of the situation, we state that (a) AA3{x,Q^) and A'E{x,Q'^) are 
relatively well known in the kinematic regions where data exist; (b) the Bjorken 
sum rule [176] is confirmed by the data; (c) the first moment of AS, and thus 
the quark-plus-antiquark spin contribution to the proton spin, is of the order 
of 25% or less (known as "spin surprise"); and (d) AAsixjQ'^) and the spin 
gluon density Ag{x,Q'^) are constrained very little by the data so far. Note 
that this finding for A^g implies also that the polarized strange density is still 
unknown to a large extent. The present situation concerning Ag is represented 



by Figure 22, which compares the polarized gluon densities of several recent NLO 
sets of spin-dependent parton distributions [^, 44, 79 1, all consistent with current 



DIS data. The wide range of possible gluon polarization expressed by the figure 
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Figure 22: The polarized gluon densities as given by six different NLO parame- 
terizations §|, ||, |7|] , at the scale Q = 10 GeV. 



does not come as a surprise. For DIS, the gluon is only determined through the 
scaling violations of the structure functions g'f'"; however, so far only fixed-target 
polarized DIS experiments have been carried out, which have a limited lever arm 
in Q^. The measurement of Ag remains one of the most interesting challenges 
for future high-energy experiments with polarized nucleons. 
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